Results for light pseudoscalars from three-flavor simulations 



The MILC Collaboration: C. Aubin? C. Bernard, a * C. DeTar, b Steven Gottlieb? E.B. Gregory? 
Urs M. Heller? J.E. Hetrickf J. Osborn, b R.L. Sugars an d D. Toussaint d 

a 

Department of Physics, Washington University, St. Louis, MO 63130, USA 
b Physics Department, University of Utah, Salt Lake City, UT 84112, USA 
c Department of Physics, Indiana University, Bloomington, IN 47405, USA 
d Department of Physics, University of Arizona, Tucson, AZ 85721, USA 
c American Physical Society, One Research Road, Box 9000, Ridge, NY 11961, USA 
f Physics Department, University of the Pacific, Stockton, CA 95211, USA 
g Department of Physics, University of California, Santa Barbara, CA 93106, USA 

We compute pseudoscalar meson masses and decay constants in partially quenched QCD with three dynamical 
flavors of improved staggered quarks. Fitting the lattice data to staggered chiral perturbation theory (SXPT) 
forms and extrapolating in quark mass and lattice spacing, we find: f w — 129.5±0.9±3.5 MeV, Jk = 156.6±1.0± 
3.6 MeV, and fic/U = 1.210(4)(13). The value for f K /U implies |V tt „| = 0.2219(26). We also obtain m u /m d = 
0.43(0)(1)(8), or, with more input on electromagnetic (EM) effects from the continuum, m u /m,d — 0.41(0)(1)(4). 
Using the one-loop perturbative mass renormalization I], we determine m u and m,d in the MS scheme. 



The properties of the tt-K system are amenable 
to high-precision lattice computation. At fixed 
lattice spacing and quark mass, our improved 
staggered |2] computation achieves a statistical 
accuracy of 0.1% to 0.4% for decay constants and 
0.1% to 0.7% for squared meson masses, down to 
m-n « 250 MeV (ra,/^ w 0.3). This allows us 
to check that the chiral behavior is as expected 
from continuum chiral perturbation theory (XPT) 
plus the controlled 0(a 2 ) corrections predicted by 
SXPT 3 . The check is especially important be- 
cause our implementation of staggered fermions 
uses the as yet unproven v^Det trick to move from 
4 "tastes" (doublers) to 1. With SXPT, the chiral 
and continuum extrapolations can then be per- 
formed to extract f v and fx, quark masses and 
ratios, and several low energy constants. A de- 
tailed account of our methods and results appears 
in Ref. Here, we give an overview and update 
our light quark mass results to take into account 
continuum analyses of EM effects. 

Our computation uses seven lattice runs; five 
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make up the "coarse" set (a 0.125 fm), with 
degenerate u, d sea quark mass am' = 0.005 to 
0.03 and strange sea quark mass am' s — 0.05, 
and two make up the "fine" set (a w 0.09 fm), 
with am' = 0.0062,0.0124, and am' s = 0.031. 
(The primes distinguish simulation masses from 
their physical counterparts.) For Goldstone me- 
son masses and decay constants we have extensive 
partially quenched data: all combinations of 8 or 
9 valence masses between 0.1m' s (0.14m' s on the 
fine set) and m' s . For other tastes, we have most 
of the full QCD pion masses but no decay con- 
stants and no partially quenched data. 

Given the available data, we can only hope to 
make precise computations for Goldstone mesons. 
However, since mesons with other tastes affect 
the Goldstone properties through chiral loops, 
we need to know the masses of those mesons to 
reasonable precision. We get the needed split- 
tings and the slope in quark mass from a tree- 
level SXPT fit to full QCD data for pions of all 
tastes. Using these tree-level masses in the one- 
loop terms gives an NNLO error in the final re- 
sults; since the maximum error in the tree-level fit 
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is only ss 7%, the error in the final results turns 
out to be well under 1%, i.e., negligible. 

To get good fits to SXPT forms, we need to 
place upper limits on the valence quark masses 
(m x , m y ) considered. We consider 3 data sub- 
sets: subset / (94 data points) with m x +m y < 
0.40m' s (coarse) and < 0.54m' s (fine); subset 77 
(240 points) with m x +m y < 0.70to(, (coarse) 
and < 0.80to^, (fine); and subset III (416 points) 
with m x +m y < 1.10to' s (coarse) and < l.\Am' s 
(fine). We can tolerate heavier valence masses 
(compared to m' s ) on the fine lattices, since both 
m' s /m s and the taste splittings are smaller. 

On subset I, the maximum valence-valence 
Goldstone mass is «350MeV. Adding on the aver- 
age taste splitting gives «500MeV; the worst case 
is «580MeV. Defining a chiral expansion param- 
eter x=(500MeV) 2 /(87r 2 /2), we thus expect the 
errors of NLO SXPT to be of order X 2 ~ 3.5%. 
Since this is much larger than our statistical er- 
rors, NNLO terms are clearly needed, even for 
this lightest mass set. Unfortunately, the NNLO 
SXPT logs have not been calculated. But for 
the higher masses where the NNLO terms are 
important, the logs should be smoothly varying 
and well approximated by NNLO analytic terms. 
Our "NNLO" fit functions therefore include the 
full NLO expressions supplemented by the an- 
alytic terms at NNLO. We fit decay constants 
and masses together, including all correlations; 
this constrains the new SXPT hairpin parameters, 
which are common to both. We also fit coarse and 
fine lattices together, which constrains the overall 
lattice spacing dependence. 

The NNLO fit has 20 unconstrained param- 
eters: 2 at LO, 4 physical ones (Li) at NLO, 
4 taste- violating ones at NLO, and 10 physical 
NNLO analytic parameters. There are an ad- 
ditional 16 tightly-constrained parameters that 
allow for the variation of physical parameters 
with a (~ a s a 2 Aq CD ps 2%). Finally, 4 tightly- 
constrained parameters allow the scale determi- 
nations to vary within statistical errors. 

We are able to find good NNLO fits for sub- 
sets I and II. In subset III, even NNLO fits break 
down. Yet it is subset III that can allow us to in- 
terpolate around m s . The strategy for subset III 
is to fix the LO and NLO terms from lower mass 
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Figure 1. Squared meson masses vs. m x /m' s . 
"Kaon" points (crosses and octagons) have m y = 
m' s , 0.8m' s , or 0.6to' s ; "pion" points (diamonds 
and squares) have m x = m y . 

fits; then add on ad hoc higher order terms to get 
a good interpolating function around m s . We use 
such fits for central values of quark masses and de- 
cay constants; results from the other subsets are 
included in systematic error estimates. Figure ^ 
shows the fit results in subset III for partially- 
quenched meson masses; while Fig. |21 shows the 
same fit for pion decay constants. Points and fit 
lines in these plots have been corrected for finite 
volume effects at one-loop in SXPT. 

To find physical values of the (bare) quark 
masses, we first extrapolate the chiral parame- 
ters to the continuum. The extrapolation is linear 
in the leading discretization errors, either a 2 s a 2 
(for taste-violating terms) or asa 2 (for "generic" 
discretization effects in the physical parameters). 
We then set valence and sea quark masses equal 
(i.e., go to "full QCD"), and iteratively interpo- 
late in m' s and extrapolate in to' until both the 
kaon and the pion reach their physical masses. 
The extrapolations of masses and decay constants 
in to' (after m' s has been adjusted to m s ) are 
shown as dashed red lines in Figs. n an d 121 

We must distinguish here between the exper- 
imental meson masses, QCD masses (no elec- 
tromagnetism) , and those with EM effects and 
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Figure 2. "Pion" decay constants with m x — m y . 

isospin violations turned off, and m#, which 
are the "physical" masses that determine rh and 
m s in our simulation. From continuum XPT, 
these masses are related by 

m\ » [(m^ D ) 2 + (m^ D ) 2 ]/2 

\ m K° > ~ m K« 

(m^ D ) 2 « m 2 Y+ -(l + A B )(m 2 + -m 2 o) , 

where m^o, m T +, m^o, and are the ex- 

perimental masses, and parametrizes the vi- 
olations of Dashen's theorem [5]. A conserva- 
tive choice, with minimal continuum input, is 
= 1 ± 1. Using more information from the 
continuum gives jS] Ae = 1.2 ± 0.5. 

Once rh and m s are known, we may continue 
the kaon extrapolation as a function of m x (green 
dashed line Fig. inset) until the value rn^ is 
reached. This determines m u up to small isospin- 
violating corrections, which arise because m u = 
and = rh for the sea quarks. 

We now extrapolate the decay constants to 
physical quark masses, obtaining 

U = 129.5 ±0.9 ±3.6 
f K = 156.6 ±1.0 ±3.8 
f K /U = 1.210(4)(13) , 



where the errors are statistical and systematic 
(coming mainly from the chiral and continuum 
extrapolations, and, for dimensionful quantities, 
the scale error). Using Ref. [7], the result for 
f K /U implies \V US \ = 0.2219(26). 

For light quark masses at 2 GeV, we obtain 

mfj =1.7(0)(1)(2)(2); 1.6(0)(1)(2)(1) MeV 

mf s =3.9(0)(1)(4)(2); 3.9(0) (1) (4) (1) MeV 

m u /m d = 0.43(0)(1)(8); 0.41(0)(1)(4) , 

where the first result in each line comes from the 
conservative choice = 1 ± 1 ; while the second 
uses Ae = 1.2 ± 0.5 ^6! . In each case the first two 
errors are from statistics and lattice systematics; 
the last, from EM effects, and the third (for MS 
masses) is an estimate of perturbative errors 

Despite the fact that we have a large number 
of free parameters in our fits, we believe that the 
chiral behavior predicted by SXPT is well tested 
by our analysis. The evidence for this is that good 
fits are not possible with any of the following al- 
ternatives: (a) fits to the continuum form (36 pa- 
rameters), (b) fits with all chiral logs and finite 
volume corrections omitted, i.e., keeping only an- 
alytic terms (38 parameters) , and (c) fits with all 
chiral logs omitted, but finite volume corrections 
included (38 parameters). In subset 77, all these 
fits have x 2 /d.o.f. > 3 with d.o.f. > 202. 
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